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Abstract 

Using scaled-particle theory for binary mixtures of two-dimensional hard particles with orien- 
tational degrees of freedom, we analyse the stability of phases with orientational order and the 
demixing phase behaviour of a variety of mixtures. Our study is focused on cases where at least 
one of the components consists of hard rectangles, or a particular case of these, hard squares. A 
pure fluid of hard rectangles has recently been shown to exhibit, aside from the usual uniaxial 
nematic phase, an additional oriented phase, called tetratic phase, possessing two directors, which 
is the analog of the biaxial or cubatic phases in three-dimensional fluids. There is evidence, based 
on computer simulation studies, that the tetratic phase might be stable with respect to phases with 
lower translational symmetry for rectangles with low aspect ratios. As hard rectangles are mixed, 
in increasing concentration, with other particles not possessing stable tetratic order by themselves, 
the tetratic phase is destabilised, via a first- or second-order phase transition, to uniaxial nematic 
or isotropic phases; for hard rectangles of low aspect ratio (hard squares in particular), tetratic 
order persists in a relatively large range of volume fractions. The order of these transitions depends 
on the particle geometry and dimensions, and also on the thermodynamic conditions of the mixture. 
The second component of the mixture has been chosen to be hard discs or disco-rectangles, the 
geometry of which is different from that of rectangles, leading to packing frustration and demix- 
ing behaviour, or simply rectangles of different aspect ratio but with the same particle area, or 
different particle area but with the same aspect ratio. These mixtures may be good candidates 
for observing thermodynamically stable tetratic phases in monolayers of hard particles. Finally, 
demixing between fluid (isotropic-tetratic or tetratic-tetratic) phases is seen to occur in mixtures 
of hard squares of different sizes when the size ratio is sufficiently large. 
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I. INTRODUCTION 



Mixtures of three-dimensional rod-like molecules have been analysed quite extensively, 
both expe—H, and t We t icaH y g BeautiM expels on — and coUoidal 
particles have been reported, with the observation of different phases, such as isotropic and 
nematic phases of various kinds {2]. Of special interest is the issue of orientational phase 
transitions and, in particular, of entropically-driven phase segregation in these systems, a 
long-debated question in the context of hard spheres (3j. There is now ample theoretical evi- 
dence for segregation and demixing phenomena in hard-core mixtures of anisotropic particles 
, |9j. Particularly interesting are segregation phenomena between two nematic 
phases, which sometimes show upper critical points jioj]. More recently, the occurrence and 
overall influence of spatially ordered phases, in the perfect-alignment approximation, have 
been analysed 0, and more complete studies, lifting the latter approximation, on the 



effect of layered phases and microsegregation phenomena in the phase behaviour have been 
carried out [12I, Q, 



Studies on the corresponding two-dimensional mixtures are very scarce. In Ref. [15] it 
was shown that two-dimensional isotropic mixtures of hard convex bodies can never demix 
within scaled-particle theory (SPT). A unimodal polydisperse mixture of hard needles was 



also studied within the Onsager approach in Ref. [161 ]. one result being that the isotropic- 
nematic transition is always of second order. Also, an equation of state for two-dimensional 
mixtures of hard bodies has been constructed starting from an approximation for their direct 



correlation functions 



171 ]; it was found that demixing never occurs. Finally, a theoretical 



study has recently been carried out for mixtures of hard rectangles and discorectangles 
within the framework of SPT 18( . Using a bifurcation analysis, demixing between different 



phases, one of which is an orientationally-ordered phase, was shown to occur. 

The analysis of this problem is sufficiently motivated by the importance of surface-phase 
transitions experienced by monolayers of adsorbed molecules. But, since the transition from 
the disordered phase to the oriented (nematic) phase in one-component fluids may be, and 
in most cases is, of second order in two dimensions, many new interesting features may 
arise in the phase diagram of the mixture, such as tricritical and critical end points, which 
are absent in three dimensions. Our recent work on two-dimensional hard-rod fluids [3] 
has demonstrated that these systems do in fact exhibit a richer phase behaviour, with the 
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additional fact that, due to spatial restrictions, phase behaviour may be more sensitive 
to subtle effects associated with the particle geometrical shapes. It would be desirable to 
understand the differences and similarities between the two dimensionalities and elucidate 
their origin. 

In fact, liquid-crystalline phase transitions in one-component fluids depend strongly on 
dimensionality. A general trend of two-dimensional systems with continuous symmetry is the 
lack of true long-range order which in a two-dimensional nematic would be reflected in 
the presence of quasi-long-range orientational order 



21 



22], 123, |24( • Thus, in the absence 



of any other mechanism, transitions from the isotropic phase (I) to the uniaxial nematic 
phase (N u ) in two-dimensional, one-component fluids of hard rods, may be governed by a 



disclination unbinding-type mechanism 



22 



23 



24j , as proposed by the Kosterlitz-Thouless 



3, 



26, 



271 ]. the nature 



(KT) theory [2_5| (in fluids of hard ellipses, for which simulations exist 
of the transition seems to depend on the aspect ratio 23] )• Mean-field theories, which cannot 
account for these effects since collective fluctuations and the dynamics of topological defects 



are not proper 
field) type [28, 



y 



or not at all) described, predict continuous transitions of the usual (mean- 
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301 ] . However, the direct role played by dimensionality in two-dimensional 



mixtures may be of secondary importance as far as entropically-driven demixing transitions 
are concerned, since local entropic effects associated with packing may completely preempt 
KT-type effects; in this sense, the mechanisms governing these systems could be intimately 
connected with those operating in the corresponding three-dimensional mixtures, which are 
described qualitatively correctly by mean-field theories of the density-functional type [lj]. 

Of particular interest is the case of hard rectangular particles, which might exhibit a 
so-called tetratic phase (N t ) in the one-component fluid |30|, |3l|, |32j, [33, [34|. The tetratic 
phase possesses two equivalent directors pointing along mutually orthogonal directions; it 
exhibits a symmetry higher than that of the particles making up the fluid. The nature of 
the transition from the I to the N t phases is unknown, but density- 



it to be of second order 



18, 



unctional studies predict 
301 ] . A recent bifurcation analysis 31], combined with further 
calculations which include three-body correlations 34|, have not been conclusive as to the 
absolute thermodynamic stability of the tetratic phase with respect to phases with spatial 
order, although preliminary computer simulations 32j, [33, [34| seem to support the tetratic 



phase as an intermediate phase between the isotropic phase and the crystalline phase, at 
least for aspect ratios less than ~ 7-9 [34|. However, questions remain as to the maximum 
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aspect ratio that can support tetratic order (for example, experiments on -non-equilibrium- 



vibrated monolayers of granular particles 35| extend it up to ~ 12), the interplay between 
the N t and the usual uniaxial nematic phase N u (possessing only one director or alternatively 
two non-equivalent directors), or the role played by the nature of the crystalline tetratic 
ohase (with an aperiodic solid being the most promising candidate 33)). Finally, recent work 



361 ] has shown the ability of these particles to promote spatial order when confined between 



parallel one-dimensional plates, with potentially interesting applications as building blocks 
for self-assembly. 

The role played by tetratic order in mixtures where one of the components consists of 
hard squares or rectangles, while the other may or may not promote this order, is the aim of 
the present paper. We investigate the phase behaviour of mixtures of hard squares or hard 
rectangles with other particles of different geometry. The ability of hard rectangles to induce 
either short- or long-range tetratic order originates from the sharp corners of their shape 
(here we do not address the question on the existence of true long-range order of nematic 
correlations in these systems). Rectangles with low aspect ratio k = L/a (with L their 
length and a their breadth) may orient along either a particular direction or the orthogonal 
direction equally easily without hampering packing efficiency, and a tetratic phase results. 
For larger aspect ratios, hard rectangles stabilise into the usual uniaxial nematic. Hard-rod 
particles terminated by a semicircle (such as a disco-rectangle) do not pack efficiently in a 
tetratic arrangement, and can only form one (uniaxial) nematic phase. When these two types 
of particles, or when two types of hard rectangles, one exhibiting a tetratic phase which is 
not present in the other, mix together, the tetratic order induced by one component survives, 
at higher pressures, in some range of particle volume (area) fractions, with an extent that 
depends on the geometrical compatibility between species. Nematic demixing occurs, and 
the resulting phase diagrams show a rich variety of features, such as phase transitions of 
different order, triple, critical, tricritical, azeotropic and critical end points. 

A particular situation is when hard squares of different particle area are mixed. A con- 
sensus is now beginning to emerge that entropy-driven demixing in additive hard-sphere 
mixtures involves spatially non-uniform phases j^]. The current situation concerning mix- 
tures of parallel (i.e. with frozen orientational degrees of freedom) hard cubes looks very 



similar 37|]. Now, despite initial evidence based on simulation that no demixing occurs in a 
mixture of parallel hard squares on a lattice [381 ]. supported by off-lattice density-functional 
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calculations on the fluid phase 391], recent simulation work gives evidence for demixing 



involving an inhomogeneous phase 42]. In contrast, our calculations, which incorporate 
orientational degrees of freedom (though, admittedtly, do not contemplate inhomogeneous 
phases) , show demixing involving two fluid phases, at least one of which is oriented (with an 
island of instability in the phase diagram and sometimes with an associated upper critical 
point), when the size ratio is sufficient large, while the completely isotropic mixture does 
not exhibit segregation. Clearly our results demand for additional computer simulations on 
the freely-rotating hard-square model, which have not been reported yet. 

After presenting a brief summary of the theoretical approach, which is a density- 
functional theory based on the SPT approximation, we present the results and conclude 
with a summary and some final remarks. The Appendix presents some details on the cal- 
culation of spinodal lines and tricritical points. 

II. THEORY 

Let us first briefly recall the theoretical model used and the approximations implemented. 
The model is based on the SPT approximation for two-dimensional binary mixtures, first 
applied in three-dimensions by Cotter and Wacker [43]. The first extension to the two- 
dimensional case, a model of hard rectangles with restricted orientations [44] . was later 
rederived for isotropic mixtures of general hard convex bodies [15]. Here we use the im- 
plementation of the SPT approximation for oriented mixtures derived, using the standard 
procedures, in Ref . [18] , to which the reader is referred for further details. Let us denote the 
free-energy density / in units of the thermal energy kT by i.e. $ = f /kT. The density 
functional for the corresponding excess (over ideal gas) quantity, $exo[^i)M' depends on 
the two orientational distribution functions h u ((f)) [y = 1,2) for the two components of the 
mixture, <fi being the angle between the particle main axis and some reference direction in 
the plane, which is arbitrarily taken as the x axis. In SPT approximation is written as [3] 

h 2 ] = p / -log (1 - V ) + «Kl 0) » I • (1) 

Here the subindeces v, r — 1,2 refer to the two components of the mixture, p = p% + p 2 is the 
total density, with p v the density of species v, 77 = p\V\ + p 2 v 2 is the total packing fraction, 
with v v the particle area for species u, and the number fractions are defined as usual by 



X V = Pv / P . AISO, Kl O) (0) = KrO) — v u — v T , where Kr(0) is the angle-dependent excluded 
volume between species v and r. The double angular average ((Kiy ) ) is defined by 



2?r /.2tt 



((v®)) = #/ (2) 



The functions ^,,-(0) are analytical; their expressions were explicitely written in Ref. [181 ]. 
except that corresponding to the cross interaction between rectangles and discorectangles 
45]. 



Note that the one-component limit of Eqn. correctly reduces to the Onsager theory 
(where only the second virial coefficient is incorporated) as the density vanishes, since SPT 
recovers the exact second virial coefficient; however, in contrast to the three-dimensional 
case, the Onsager theory is not rigorously correct in the hard-needle limit in two dimen- 
sions, since three- and higher-order virial coefficients are not vanishingly small in this limit, 



some of them being even negative 46(). Therefore, higher-order contributions in density 
are only approximately reproduced by Eqn. ([1]). SPT can be considered as a sophisticated 
Onsager theory, in that spatial (not orientational, which are still included up to second 
order in density) correlations are somehow resumed into density-dependent terms. Al- 
ternative approaches, such as a two-dimensional version of the Parsons-Lee approach for 
three-dimensional hard rods, have the same structure and, for lack of a detailed performance 
analysis, can in principle be considered to be equivalent. The ideal contribution 

Qu&iM) = X)^^logp„-l + log [2irh v (4i)]j , (3) 

is added to obtain the complete free energy-functional hz] = $id[^i> ^2] + ^oxj^i, ^2]- 
Functional minimisation of /12] with respect to h u (4>) gives the equilibrium configuration 
of the mixture. The pressure follows from the equation 

' UT 



Instead of numerically solving the Euler-Lagrange equations associated with <3>[/ii, /12], 
we follow common practice in our research group and tackle the direct minimisation of the 



functional. Various strategies are possible 43]. In the present article we choose to intro- 



duce a parameterized form for the orientational distribution functions, h v {<p). To facilitate 
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TABLE I: Possible values of the variational parameters A^ and order parameters in the 
isotropic (I), uniaxial nematic (N u ) and tetratic (N t ) phases. 



computations, these functions are simply parameterised as 

p K ( J- ] cos 2</>+AS 2> cos 4$ 

hM = 



2- 



//gA^ 1 ' cos 2^'+Ai 2) cos 4(f>' 



(5) 



The parameters A£ take care of the two types of orientational symmetries, either uniaxial 



(k = 1) or tetratic (k = 2). Equivalently two order parameters, qi h \ can be defined as 



2rr 



d(ph u (<j)) cos (2k(j)), k = 1,2, 



(6) 



which are proportional to the coefficients of a Fourier-expansion of the functions h v (<p) 
including the two lowest symmetries that a rectangular particle can generate. Table [I] sum- 

(k) (k) 

marises the different phases with their associated values of the A u and ql parameters 
in a one-component phase. The equilibrium configurations of the mixtures are more con- 
veniently obtained by minimising the Gibbs free energy per particle g = (p + f)/p with 
respect to the parameters at fixed value of the pressure p and composition x = X\ (we 
will henceforth arbitrarily associate x with the number fraction of the component labelled as 
1). Minimisations were performed using an efficient routine based on the Newton-Raphson 
method. Coexistence (binodal) lines were located by means of a standard common-tangent 
construction. 

The use of the parameterisation ([5]) obviously introduces an approximation over the exact 
calculation, and it would be of interest to know the amount of error introduced. As far as 
the calculation of bifurcation or spinodal lines of the various phase transitions is concerned, 



the parameterisation has no impact, since first-order terms in A 



(fc) 

V i 



hu{<P) = 7T- (1 + Al 1} cos 20 + A[, 2) cos 40) + O (A u f 

Z7T 



(7) 
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which appear as quadratic terms in the free energy, are treated exactly. However, binodal 
lines and tricritical points are affected by the parameterisation, since their location depend 
on higher-order terms in A^. We have checked the parameterisation in two ways: first, some 
selected calculations have been performed using an additional cosine term, cos 60 (cos 80), 
in the parameterisation for the uniaxial (tetratic) nematic phase; the differences found, at 
the level of coexistence lines, were found immaterial. Second, tricritical points have been 
evaluated exactly by computing the exact fourth-order terms (involving the above cosine 
terms depending on the symmetry of the phase). In all cases the differences with respect 
to the calculations using (JHJ) have been found to be of minor importance. The analysis is 
presented in the Appendix, which contains quantitative details on this issue. 



III. RESULTS 



Results are presented, in the form of pressure-composition phase diagrams, in Figs. [T]- 
IH1 We divide the presentation by first showing results for different mixtures containing 
squares and discs, followed by mixtures of hard rods (rectangles and discorectangles). In the 
following we will generally use the following acronyms for the different particles: HS (hard 
squares), HD (hard discs), HR (hard rectangles) and HDR (hard discorectangles). 



A. Mixtures of hard squares 

The fluid of freely-rotating hard squares has been studied by Monte Carlo simulations 
321 ]. with indications that a tetratic phase might be stable prior to crystallisation. On 



the other hand, mixtures of hard squares have been analysed by computer simulations 



38, 



421 ] and density-functional theory [39J] in the approximation of perfect orientational order 



Indications that there is no demixing in this system 27], |38( were later challenged by evidence 



for a spinodal line from Monte Carlo simulation [42| ; the demixing transition involves a fluid 
phase and a phase with big squares in close-packed aggregates. Since orientational disorder 
may be important for these mixtures, we believe it is of interest to address this problem with 
the SPT approximation, where orientational disorder is allowed (obviously our approach can 
only search for demixing behaviour involving fluid, either isotropic or tetratic, phases). 
The results shown in Fig. [1] correspond to a mixture of freely-rotating squares with size 
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FIG. 1: Phase diagram of a hard-square (L„ = a v for v = 1,2) mixture with K4 = 10 and K2 = 1 
in the scaled-pressure pv2/kT vs. composition x = xi plane. Open circles: critical points. Filled 
circles: critical end points. Inset is a qualitative scheme depicting the relative locations of the I-N t 
spinodal and the region of demixing. Two-phase region is indicated by the gray area. Note that 
the scale of the composition (horizontal) axis is discontinuous. 

ratio 1:10, and have been gathered in a pressure-composition phase diagram, with label 
1 assigned to the larger particles. According to SPT theory, the fluids of the unmixed 
(one-component) species both exhibit corresponding second-order I-N t phase transitions 
at reduced pressure pv/kT ~ 0.52 (here v is the particle proper area). Note that, for 
squares, a uniaxial nematic phase is not possible by construction. For the mixture, the 
only possible (fluid) phases are also the I and N t phases. Remarkably, the mixture exhibits 
segregation between two tetratic phases with different area fractions occupied by the two 
components. The region of demixing is a closed loop bounded above (below) by an upper 
(lower) critical point (in fact, since the I-N t spinodal line crosses the demixing region -see 
inset showing schematic topology-, segregation mostly proceeds between I and N t phases). 
Since mixtures of perfectly parallel hard squares do not demix it is surprising that, in this 
mixture, orientational disorder induces demixing. 

Closed loops of immiscibility have been predicted in mixtures of hard rods in the Onsager 
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approximation [7|, [8| (see however Ref. [481]). and also in fluids of parallel hard rods [49] , 



both in three dimensions. Our calculations show that they are also a property of some two- 
dimensional mixtures. In our case demixing occurs at fairly low values of composition; this 
is easy to explain since we expect segregation to take place when the volume fraction of both 
components are approximately equal, rji ~ r^ 2 , which implies that segregation will occur for 
x ~ V2/V1 = 10~ 2 . As the size ratio of the squares decreases, the island of immiscibility 
diminishes, and we are eventually left with a second-order I— N t transition in the whole 
composition interval. The demixing island disappears when the size ratio is equal to 1 : 4 
(see Appendix). 

The mechanism underlying N t -N t demixing is different for the upper and lower parts of 
the demixing island. As is well known, demixing and ordering phenomena in mixtures of 
hard anisotropic particles result from the competition of entropies of different origin: mix- 
ing, orientational and excluded-volume entropies. In the transition from the isotropic to the 
nematic phase the last two terms compete. In mixtures, the excluded-volume entropy has 
contributions from the two species and from the unlike-particle interactions. The mechanism 
explaining the lower N t -N t segregation phenomena in Fig. [I] is the classical one: the balance 
between mixing entropy and the excluded-volume term of unlike species, which counter- 
balances the tendency toward mixing of the former. However, the upper N t — N t segregation 
region has a different origin, since here orientational entropy plays a role. This might have 
been suspected a priori, since in the limit of perfect order (parallel squares) no demixing 
occurs. At sufficiently high pressures in the freely-rotating fluid, sufficiently near the close- 
packed limit, the orientational order is almost saturated and can no longer compete with the 
other terms, so that no demixing is expected. As pressure is reduced, orientational entropy 
begins to play a role and in fact the contribution from the small squares is the driving force 
toward demixing. 

For the size ratio shown (1:10), two critical end points, indicated in the figure by filled 
circles, appear in the phase diagram, defined by the points where the I-N t spinodal line 
crosses the demixing region; the location of these points, given by x^* = 0.012, p( 1 '*V2/kT = 
13.86 and x^* = 8 x 10 -5 , p( 2 )*t> 2 /fcT = 49.38, can be approximated from the set of 
candidates to tricritical points obtained by means of a bifurcation analysis (here and in what 
follows, numerical values for the location of tricritical and points should be understood to 
result from rigorous bifurcation analysis of the free-energy functional, and may in some cases 
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FIG. 2: Phase diagram for a HS/HD mixture in the scaled-pressure pv\/kT vs. composition x = x\ 
plane, with v± = L\U\ the volume of the squares. The side length of the squares, L\ = a\, is chosen 
to be the same as the diameter of the discs, oi- Filled circle indicates tricritical point. Two-phase 
region is indicated by the gray area. 

be at variance with those obtained from the variational minimisation, on which all the phase 
diagrams presented are based; see Appendix for details on the bifurcation analysis). As the 
size ratio is diminished, first the upper critical point disappears, and the upper critical end 
point becomes a tricritical point; then the lower critical end point becomes a tricritical point, 
and the lower critical point also disappears, leaving two tricritical points before the whole 
demixing region vanishes. The lower critical point always stays within the N t region (i.e. no 
II demixing is observed), . 

Demixing in this mixture is therefore associated with a symmetry-breaking phase, the 



tetratic phase. There are strong arguments 15j disallowing demixing in the I phase (within 



the context of the SPT approach). This would not necessarily imply that demixing is com- 
pletely ruled out in the rotationally-symmetric I phase of hard-square mixtures or mixtures 
of particles with different geometries, but a more sophisticated theory, incorporating exact 
higher-order virial coefficients (which are known to be important for two-dimensional hard 
convex bodies), would be necessary to settle this point. 
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B. Mixtures containing squares and discs 

Fig. [2] shows the phase diagram corresponding to a mixture of hard discs and squares 
(the diameter of the former being equal to the side length of the latter). Here only one 
of the components (the squares) has tetratic order. By choosing the particle areas of both 
species to be approximately equal (the ratio being ~ 0.79) we focus on the role of parti- 
cle geometry. In this mixture we find that the N t phase exhibited by the pure system of 
squares survives up to a maximum concentration of discs of ~ 50% (this corresponds to an 
area ratio approximately equal to 0.79). Complete demixing occurs at high pressure, while 
the I— N t transition becomes of second order below a tricritical point, located at x = 0.57, 
pvi/kT = 131.95 (see Appendix). Fractionation becomes stronger as pressure is increased, 
since excluded- volume considerations are increasingly important in this limit: geometri- 
cal mismatch (given by the unlike disc-square interaction) grossly counterbalances mixing 
entropy, inducing strong segregation. 



C. Mixtures of rectangles and discorectangles 



Hard rectangles have been predicted to exhibit a phase with tetratic order when their 



aspect ratio is sufficiently low 



30 



31 



32 



33 



34( | . Tetratic order is enhanced as the aspect 



ratio of the rectangles is decreased. However there is some uncertainty as to the critical 
value of k beyond which the tetratic phase is no longer possible; the SPT approach indicates 
that for k < 2.21 the stable nematic phase is the tetratic nematic 50j (as mentioned in 



the introduction, computer simulations indicate that this value may be much larger, in the 
range ~ 7-9). In contrast, a fluid of hard discorectangles can only support a uniaxial nematic 



phase [23 ]. 



In previous work 18j we analysed possible demixing scenarios of mixtures of hard rect- 
angles (HR) and mixtures of hard discorectangles (HDR), using the same SPT theory. In 
this section we further investigate this problem by considering a wider range of mixtures, in 
particular, crossed mixtures (i.e. mixtures consisting of HR and HDR particles); the effect of 
particle geometry is an aspect that can be assessed in a very direct way by analysing crossed 
mixtures, as well as the role played by the tetratic phase against the standard uniaxial 
nematic phase when particles of different geometries are mixed. Results will be presented 
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X 

FIG. 3: Phase diagram for a HR/HDR mixture in the scaled-pressure pv\/kT vs. composition 
x plane. The two components have the same aspect ratio {n\ = L\ju\ = 2 for the rectangles 
and k-2 = (L2 + <?2)/&2 = 2 for the discorectangles), and same particle area. Open circle indicates 
critical point, shaded circle denotes an azeotropic point, while filled circle indicates tricritical point. 
Two-phase regions are indicated by gray areas. 

by means of phase diagrams including, not only spinodal lines, but also binodal lines when 
present. 

In the following we first consider mixtures of hard rods, one of which can be stabilised 
into a tetratic phase (i.e. HR particles with aspect ratio k < 2.21) while the other cannot 
(i.e. either HR particles with k > 2.21, or HDR particles). We intend to understand how 
the tetratic phase is destabilised by the geometrical mismatch of the particles. 

The first mixture that we consider is a mixture of HR and HDR particles, both with 
k = 2 [in the latter case the aspect ratio parameter is defined as k — {L + a) j 'a]. Also, the 
same particle areas have been chosen for both types of particles, in an attempt to single out 
features of phase behaviour mainly driven by differences in particle geometry. The resulting 
phase diagram is depicted in Fig. [3j At high pressure a demixing region occurs between two 
uniaxial nematic phases (i.e. there is N u -N u demixing), bounded by a lower critical point. 
Fractionation becomes stronger as pressure increases, an effect ultimately associated with the 
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FIG. 4: Phase diagram for a HR/HDR mixture in the scaled-pressure pv\/kT vs. composition 
x plane. Values of the parameters are k\ = 1.5, a\ = 1 for the rectangles and K2 = 2 and same 
particle rectangle of aspect ratio equal to 2 and unit breadth for the discorectangles. 

Open circle indicates critical point, shaded circle denotes an azeotropic point, while filled circles 
indicate tricritical points. Two-phase regions are indicated by the gray areas. 

slight difference in particle geometry, which penalises the mixed state due to unfavourable 
packing between dissimilar particles; this arises from the circular caps of the HDR particles. 
At lower pressure there is a transition between the isotropic and the uniaxial nematic phases, 
which occurs via a first-order phase transition with small fractionation. The phase diagram 
is of the azeotropic type: an azeotropic point (no fractionation) appears at x ~ 0.32, with an 
associated (small) pressure range where the N u phase is reentrant. Since the HR component 
of the mixture possesses a stable tetratic phase, an island of tetratic order exists for high 
values of x, separated from the isotropic phase via a second-order phase transition and from 
the uniaxial nematic via a first-order transition (of course this is also the case for the pure 
HR fluid, x = 1). The island of N t stability is very small; for a HR fluid with k = 2, the 
range of pressures where the tetratic phase is stable is small |31[ but, in addition, tetratic 
order is easily destroyed when adding to the mixture particles that do not conform with 
tetratic symmetry. 
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FIG. 5: Phase diagram for a HR/HDR mixture in the scaled-pressure pv\/kT vs. composition 
x plane. Values of the parameters are K\ = 2, a\ = 1 for the rectangles and k 2 = 1.5 and same 
particle area as rectangle of aspect ratio equal to 1.5 and unit breadth for the discorectangles. 
Filled circle indicates critical end points. Two-phase region is indicated by the gray area. 

Next we consider two different mixtures that represent slight variations with respect to 
the previous mixture. In Fig. H]we analyse a HR/HDR mixture, where the HR component 
has been shortened to an aspect ratio of K\ = 1.5, keeping the breadth to the same value, 
while the HDR component still has k 2 = 2 but the particle area is made equal to that of 
the HR particle of the previous mixture. Essentially we would like to analyse the effect 
of shortening the hard-rectangle component in the previous mixture (while at the same 
time reducing its particle area by 25%) with the aim of increasing the strength of tetratic 
ordering. Since the tetratic phase of the one-component HR fluid is much more stable, the 
tetratic phase in the mixture stabilises into a larger range of compositions. The diagram 
is topologically equivalent to that in Fig. [31 except that the critical end point in the I-N t 
spinodal line now becomes a tricritical point, and there appears a (new) I— N u — N t triple point. 

Non-trivial changes are obtained if the mixture of Fig. [3] is changed by shortening the 
HDR particle down to an aspect ratio of k% = 1.5 while the particle area is made equal to 
that of a HR particle of aspect ratio 1.5 and unit breadth; as a result, the nematic phase 
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X 



FIG. 6: Phase diagram for a HR/HR mixture in the scaled-pressure pv\/kT vs. composition x 
plane. Particle breadths are set to the same value, while aspect ratios are K\ = 10 and k<i = 5. 
Open circle: critical point. Filled circle: tricritical point. Two-phase regions are indicated by the 
gray areas. 

of the one-component HDR fluid appears at a much higher pressure. The phase diagram is 
presented in Fig. [H This phase diagram can be regarded as a continuation of that in Fig. 
[3] where the N u -N u demixing region and the I-N u coexistence have collapsed into a single 
demixing region. The I-N t phase stability region now increases, as a result of the HDR 
particle having a lower proper area. 

D. Mixtures of rectangles 

Finally we consider mixtures of HR particles. Here the breadth of all particles will be 
taken to be unity, and we change the length or, equivalently, the aspect ratio. Three cases are 
considered: (i) K\ — 5, k% — 10, so that none of the components exhibits tetratic symmetry; 
(ii) K\ = 2, k 2 = 1-5, with both species having tetratic phases; and (iii) K\ in the range 
4.0-5.0, k,2 = 2, so that only the second species can stabilise into a tetratic phase. This 
cases are shown in Figs. EHHl respectively. 

In mixture (i) no tetratic phase appears (Fig. [6]). There is N u -N u demixing at high 
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FIG. 7: Phase diagram for a HR/HR mixture in the scaled-pressure pv\/kT vs. composition x 
plane. Particles have the same breadth, with k\ = 2 and k<i = 1.5. Two-phase region is indicated 
by the gray area. 

pressure. The I-N u transition is of first order, but becomes of second order at a tricritical 
point (see inset), located at x — 0.35, pvi/kT = 4.07 (see Appendix). 

The aspect ratios of the two components of mixture (ii), K\ — 2 and Ki = 1.5, were chosen 
such that both possess a stable tetratic phase. This requires their values to be very similar; 
therefore, no demixing occurs (Fig. [7]), and a relatively featureless phase diagram results. 

Finally, in mixture (iii), the shorter component exhibits a tetratic phase that propagates 
from the x = axis into finite values of composition when the longer component is added; 
however, due to the considerable difference in lengths between the two species, the N t phase 
does not survive very much as a function of composition; in fact, the region of tetratic 
stability can hardly be seen in the phase diagram, Fig. [HI At high pressure there is a large 
demixing region which, as the particle aspect ratio of the first component is increased (from 
4.0 to 5.0), expands considerably. Again demixing is due to unfavourable excluded-volume 
interactions between unlike species. In general, the phase diagram varies considerably even 
for very slight variations in particle shape of the first species. For n\ = 4.0 and 4.6, a 
separate region associated with the first-order I-N u transition appears at low pressures; in 
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FIG. 8: Phase diagram for a HR/HR mixture in the scaled-pressure pt^/fcT vs. composition x 
plane. All rectangles are taken to have the same breadth, and K2 = 2. The first component has 
K\ = 4.0, 4.6 and 5.0, corresponding to two-phase regions in lighter gray in the graph. 



some cases (cf. the case K\ = 4.6) reentrant behaviour in the nematic phase (stronger than 
in three-dimensional hard-rod mixtures -see e.g. Ref. is found. Also, for the case 

Ki = 4.6 an upper critical point associated with N u -N u demixing can be seen. This feature 
las been seen in three-dimensional mixtures of thin and thick hard rods of the same length 



5l| ; in our two-dimensional mixture the particle thickness is the same, and it is particle 



length that is different. Therefore, in some range of parameters two N u -N u coexistence 
regions can be found, with upper and lower critical points. These two regions merge as K\ 
increases, giving rise to a very large demixing region. 



IV. DISCUSSION AND CONCLUDING REMARKS 

To conclude, we have analysed the phase behaviour of mixtures of two-dimensional hard 
rods, using the scaled-particle theory approximation. We have particularly concentrated 
on mixtures containing hard rectangles as one of the components, with the other being a 
hard discorect angle, what we have called crossed mixtures. The geometry of hard rectangles 
allows for a tetratic phase to be stabilised; addition of a second component which does not 
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support this symmetry obviously tends to destroy tetratic order but, when the particle areas 
of both components are not very dissimilar and the pure fluid of hard rectangles strongly 
stabilises the tetratic phase (i.e. when their aspect ratio is relatively low), tetratic order 
may survive even for large area fractions of the non-tetratic-forming component. The case 
of a mixture of hard rectangles, one with tetratic order and the other without, shows the 
same trends. 

As a special case we have also considered mixtures of hard squares, and found that these 
mixtures do exhibit nematic demixing, contrary to the case where rotational degrees of 
freedom are frozen; in the latter case the mixture never segregates into two fluid phases 
(though simulations seem to point to segregation between fluid and non-uniform phases), 
while in the former there is a critical size ratio (1:4) above which there is segregation. 
This is a remarkable case where orientational disorder induces entropy-driven demixing (i.e. 
mixing order). Inclusion of non-uniform phases in our theoretical analysis would certainly 
be interesting and is left for future work. Growing evidence indicates that demixing in 
additive hard-core mixtures always involves at least one inhomogenous phase. An interesting 
possibility is that freely-rotating hard-core mixtures (which can be viewed as a special case 
of non-additive mixture) may in some cases segregate into two fluid phases, provided nematic 
phases are involved; it could be that some type of symmetry breaking, either positional or 
orientational, is required for demixing to occur. 

In general, phase diagrams of two-dimensional hard-rod particles closely resemble those 
of three-dimensional mixtures, save the fact that the isotropic-nematic transition in two 
dimensions may be (and actually is in most cases) of second order, which adds an element of 
complexity to the corresponding phase diagrams. Nematic-nematic demixing generally oc- 
curs in these systems and, for sufficiently large size ratios between the components, demixing 
competes with the isotropic-nematic (either uniaxial or tetratic) transition. Entropic effects 
due to balance between excluded-volume interactions of like and unlike components and 
mixing entropy produce strong fractionation. These effects probably counteract collective 
fluctuations leading to KT-type behaviour; for example, within our mean-field density- 
functional treatment, the continuous isotropic-nematic transition of the one-component 
fluid generally continues as a critical transition in the mixture, but sooner or later, beyond 
some value of particle composition, a demixing (first-order) transition is met via a tricritical 
point or otherwise (critical end point). However, the true nature of the isotropic-nematic 
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transition in these mixtures, and of the uniaxial and tetratic nematic phases themselves (i.e. 
whether they possess true or quasi-long range order) is an issue that should be investigated 
by means of detailed computer simulations and possibly also by experiments on vibrated 
monolayers. 

In this work we have not considered spatially-ordered phases, such as smectic, columnar 
or crystalline. Certainly some of these phases will appear at some pressure, and probably 
some of the phase behaviour shown for the mixtures considered will be preempted by these 
phases. More sophisticated density-functional treatments are required to assess this point; 



some 



terns 



proposals have already been done for one-component two-dimensional hard-rod sys- 



31] , but even the consequences of these approaches have not been explored yet; this 



will be the subject of future work 



52]. 
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APPENDIX A: BIFURCATION ANALYSIS 



In Ref. 18J two of us carried out a bifurcation analysis to study the nature of phase 
transitions in two-dimensional hard-rod fluid mixtures and to calculate the location of the 
tricritical points present in their phase diagrams. We refer the reader to this work for details 
on the calculations. For the sake of completeness, a brief summary of the main ingredients 
of the bifurcation analysis is presented here. 

The Fourier series representation of the orientational distribution functions of the two 
different species, labelled as v— 1,2, is 



h, 



1 
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h t ] cos(2fc0) 



k>l 



(Al) 



After inserting these expressions in the free-energy per particle = $/p and expanding in 
Taylor series with respect to the Fourier amplitudes h% > we obtain a Landau expansion for 
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Atp = ifiN — fi, the free-energy difference between the orientationally ordered and isotropic 
phases, in terms of these amplitudes. Further, minimizing Aip ({h^} \ with respect to all 
amplitudes except one (here chosen as hi = hf\ with i = 1,2 for uniaxial and tetratic 
nematic phases, respectively), and substituting the result back in A<p, keeping terms up to 
fourth order, we obtain the expression 

= Ah 2 + Bht, (A2) 

The coefficients A(x,rj) and B(x,rj) are both functions of the composition x = X\ and 
packing fraction 77. The spinodal curve of the transition between the isotropic and nematic 
(uniaxial or tetratic) phases can be calculated as A(x,rj*) = 0, which defines the packing 



unction of the composition. 



fraction 77* (x) as a 

It can be shown 18j that the instability region of the mixture with respect to composition 
and volume fluctuations defines a region in the plane x), at fixed K2, bounded by curves 
calculated as the roots of the function Tjy, defined by 
1 



where 



2B* 



(A3) 



T, = - (y 2 ^-) f 9 ^ 1 ) - ( 92(fI \ ( \ I ) 

dy \ dy J \ dx 2 J \ dxdy J 

Here y = p/(l — rj), and the asterisk over any function of y means that this function is to 

be evaluated at y*, its bifurcation value (note the different definition of y with respect to 



the one used in Ref. 18] ; both coincide when v u — 1, the constraint used in Ref. 18j). 
Finally, A* and A* are the partial derivatives of A with respect to x and y, respectively, 
both evaluated at y*. 



1. Mixtures of hard squares 



The packing fraction at the spinodal line obtained as the solution of A(x, r)*) = is a 
constant, independent of x, and is equal to if = [1 + 8/ (157r)] _1 . Further, the 1-N t tricritical 
points can be calculated as the roots of 



T, 
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xix 2 



1 - TT S ^ 
11 - 



1 + s 2 2 + 15^- 



'1 J 



(A5) 
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with 




1,2, 



S3 



16 (a 4 ) (a 2 ) -5 (a 3 ) 2 ' 



(A6) 



where we defined (u 11 ) = Yli x i u t- Introducing the new variable £ = X2V2/(v), and after 
some lengthy but straightforward calculations, we arrive at 

16 (25^ 2 + (34 - 16r)z + 9) 



(A7) 



x\X2 (11 + z(16r + 6 — 5z)) 
where the new variable z = (r — 1)£ was introduced, and we defined r = o\joi. The 
quadratic equation of z given in the numerator of ( 1A7I) has the following roots: 

6,2 = 25(r X _ 1} [8r - 17 ± V(r-4)(4r-l)] . (A8) 
We can see from this expression that demixing can only occur for r > 4. The pressure at 



the tricritical point can be calculated using the value of the roots found above, resulting in 



kT 



15vr 



m + y [r - (r - m 



For r = 10 we find 



£ 12 = ^(21±4V26 



(A9) 



(A10) 



with the approximate values £i ~ 0.5519 and £2 — 0.0081. Using x^ = £/ [r 2 — (r 2 — 1)£], 



we find x 



(i)* 



0.0122 and x 



(2)* 



8.12 x 10 , while the values for the pressures are 



P 



al/kT ~ 13.8603 and p^aj/kT ~ 49.3842. 



In Fig. [9] are shown the mixed and demixed states of the hard-square mixture in the 
£-r -1 plane, as obtained from the roots in Eqn. ( 1A8I) . We have checked that this is indeed 
the scenario for values of r not too much larger than its critical value r* = 4, i.e. the 1-N t 
transition changes from second (first) to first (second) order above (below) the lower (upper) 
tricritical point. As we have already shown, the phase diagram for r = 10 has two critical 
points, but the tricritical points are located close to the critical end points (see Fig. [[]). 



2. Mixtures of hard squares and hard discs 



The solution to the equation A(x, rj*) = gives us the following expression for the 
isotropic-tetratic spinodal curve: 



V 



1 + 



8xaf ' 
15tt(v) 



(All) 
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FIG. 9: The mixed and demixed states that follow from the solution to Tjy = 0, with T^ t given 



■N t 



by Eqn. flXg) - 



while the solution to TjjV = can be found also analytically as 



x 



;| 16 f Affx 



35 \n(J2 



(A12) 



Finally, the pressure at the tricritical point can be calculated as 



kT 



15tt 



5X 



1 + 



15 



2x*af 



7roy 

(Ji — | X" 



(A13) 



L 4 V 4 

For G\ = (?2 — 1 we find [x* ,p*v\/kT) ~ (0.5744, 131.9459). The value for the composition 
calculated from the minimization is about 0.53 (see Fig. [2]). This difference is due to 
the parametrization used. One should include, in the exponential parameterization of the 
orientational distribution functions, terms such as cos 40 + cos 80 to properly take 
into account the tetratic symmetry about the tricritical point. 



3. Mixtures of hard rectangles 



For mixtures of hard rectangles the isotropic-uniaxial nematic spinodal calculated as the 
solution to A(x, r)*) = gives 



1 + 



2 {{L-afY 
3n (v) 



(A14) 
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The value of the pressure at bifurcation is 



p v 1 



y 



1 + 



3((L + a))< 



(A15) 



2{(L-af)_ 

Solving Tn(x,i]*) = numerically with respect to x for a particular mixture with (Li = 
10, <7i = 1) and (L 2 = 5, er 2 = 1), we find a value of ~ 0.3472, which defines the location 
of the tricritical point. The packing fraction and the pressure at this point are rf ~ 0.4515 
and p*v\/kT ~ 4.0659, respectively. 
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